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For closed quantum systems driven away from equilibrium, work is often defined in terms of pro¬ 
jective measurements of initial and final energies. This definition leads to statistical distributions of 
work that satisfy nonequilibrium work and fluctuation relations. While this two-point measurement 
definition of quantum work can be justified heuristically by appeal to the first law of thermody¬ 
namics, its relationship to the classical definition of work has not been carefully examined. In this 
paper we employ semiclassical methods, combined with numerical simulations of a driven quartic 
oscillator, to study the correspondence between classical and quantal definitions of work in systems 
with one degree of freedom. We find that a semiclassical work distribution, built from classical 
trajectories that connect the initial and final energies, provides an excellent approximation to the 
quantum work distribution when the trajectories are assigned suitable phases and are allowed to 
interfere. Neglecting the interferences between trajectories reduces the distribution to that of the 
corresponding classical process. Hence, in the semiclassical limit, the quantum work distribution 
converges to the classical distribution, decorated by a quantum interference pattern. We also derive 
the form of the quantum work distribution at the boundary between classically allowed and forbid¬ 
den regions, where this distribution tunnels into the forbidden region. Our results clarify how the 
correspondence principle applies in the context of quantum and classical work distributions, and 
contribute to the understanding of work and nonequilibrium work relations in the quantum regime. 
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I. INTRODUCTION 

Work is a familiar concept in elementary mechanics and a central one in thermodynamics. In recent years, interest 
in the nonequilibrium thermodynamics of small systems [1-4] has motivated careful examinations of how to define 
quantum work [5-47]. In this context one often considers a process in which a quantum system evolves under the 
Schrodinger equation as its Hamiltonian is varied in time - for instance a quantum particle in a piston undergoing 
compression or expansion [48]. It is typically assumed that the system is initialized in thermal equilibrium, and the 
question becomes: how do we appropriately define the work performed on the system during a single realization of 
this process? 

One answer involves two projective measurements of the system’s energy, at the start of the process (t = 0) and 
at the end {t = r). If the system Hamiltonian is varied from H{0) = Ha to H{t) = Hb, then the measurement 
outcomes will be eigenvalues of these operators. The work performed during the process is then defined to be the 
difference between these two values, e.g. 

W = E^-Ei ( 1 ) 

if the measurements produce the m’th and n’th eigenvalues of the initial and final Hamiltonians. In this definition, work 
is inherently stochastic, with two sources of randomness: the statistical randomness associated with sampling an initial 
energy E^ from the canonical equilibrium distribution (Eq. 21), and the purely quantal randomness associated with 
the “collapse” of the final wavefunction upon making a projective measurement of the final energy (Eq. 22) [49]. 

In an analogous classical process, a system is prepared in thermal equilibrium, then it evolves under Hamilton’s 
equations as the Hamiltonian function is varied from Ha{z) to Hsiz), where z denotes a point in the system’s phase 
space. Since the system is thermally isolated, it is natural to define the work as the change in its internal energy: 

W = Heizr) - Ha{zo) (2) 

for a realization during which the system evolves from zq to Zt- As in the quantum case IE is a stochastic variable, 
but here the randomness arises solely from the sampling of the initial microstate zq from an equilibrium distribution. 

Much of the recent interest in quantum work has been stimulated by the discovery and experimental verification of 
classical nonequilibrium work relations [50-55], which have provided insights into the second law of thermodynamics, 
particularly as it applies to small systems where fluctuations are important [2]. For thermally isolated systems, 
the quantal counterparts of these relations follow directly from the definition of work given by Eq. 1 [6, 7, 9, 15]. 
Moreover, experimental tests of quantum nonequilibrium work relations have recently been proposed [56-59] and 
implemented [60, 61], providing direct verification of the validity of these results. 

Despite the evident similarity between Eqs. 1 and 2, and despite the relevance of Eq. 1 to nonequilibrium work 
relations, the definition of quantum work given by Eq. 1 might seem ad hoc. Indeed, given the absence of a broadly 
agreed-upon “textbook” definition of quantum work, one might reasonably suspect that Eq. 1 has been introduced 
precisely because it leads to quantum nonequilibrium work relations. Furthermore, the wavefunction collapse that 
occurs upon measuring the final energy has no classical counterpart. Due to this collapse, the quantum work has 
no independent physical reality until the measurement is performed at t = r. By contrast the classical work can be 
viewed as a well-defined function of time, namely the net change in the value of the time-dependent Hamiltonian 
along the trajectory Zt- Our aim in this paper is to use the tools of semiclassical mechanics, together with numerical 
simulations, to investigate the relationship between Eqs. 1 and 2, and specifically to clarify how the correspondence 
principle applies in this context. 

We restrict ourselves to systems with one degree of freedom, for which there exist explicit semiclassical approxi¬ 
mations of energy eigenstates. We focus on the transition probability P^{n\m) from the m’th eigenstate of Ha to 
the n’th eigenstate of Hb, and on its classical analogue, P'"{n\m), defined in Sec. H. In Sec. HI we investigate both 
quantities numerically for the example of a forced quartic oscillator. For high-lying initial energies P'^ oscillates 
rapidly with n, whereas P^ is a smooth function over a finite range nmin < n < Umax (Fig. 2). After integrating 
over the oscillations the two quantities are nearly identical (Fig. 3), which provides some justification for viewing 
Eq. 1 as the quantal counterpart of Eq. 2. In Sec. IV we derive P^^{n\m), a semiclassical approximation for P^ 
expressed as a sum over classical trajectories. Each trajectory in this sum carries a quantal phase, giving rise to 
coherent interferences between trajectories. When these interferences are neglected, P^'^ agrees with when they 
are included, P^'^ accurately captures the oscillations in P^ (Fig. 5). Thus the oscillations in P‘5 can be understood 
as a quantum interference pattern superimposed on a classical background. This picture breaks down around rimin 
and Umax, where P‘5 exhibits tails that tunnel into classically forbidden regions. In Sec. V we derive a semiclassical 
approximation, expressed in terms of the Airy function, that accurately describes these tails (Fig. 9). 

The analysis in Secs. IV and V relies on theoretical tools that are familiar in the field of semiclassical mechanics. 
To the best of our knowledge, these tools have not previously been applied to study the relationship between classical 
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and quantum work distributions, although similar analyses have been performed in the context of molecular scattering 
theory [62-64] and laser-pulsed atoms [65, 66], among other examples. In particular, our calculations in Sec. IV B 
closely parallel those of Schwieters and Delos [66]. 


II. CLASSICAL AND QUANTUM TRANSITION PROBABILITIES 

Here we introduce notation and specify the problem we plan to study. In Sec. IIA we define a discretized classical 
work distribution, Eq. 17, that can be compared directly to the quantum work distribution, Eq. 20. These distributions 
involve classical and quantum transition probabilities, P'^{n\m) and P'^{n\m)^ which will be the central objects of 
study throughout the rest of the paper. 


A. Classical setup 

Consider a system with one degree of freedom, described by a Hamiltonian 

H{z-X) = ^ + V{q-X) , (3) 

where z = (g,p) denotes a point in phase space, and A is an externally controlled parameter. We assume that the 
energy shells (level surfaces) of the Hamiltonian form simple, closed curves in phase space. We will consider the 
evolution of this system under the time-dependent Hamiltonian H{z] At), where A is varied from Aq = A to A,- = B. 
For compact notation, we define Ha{z) = H{z;Xo) and Hb{z) = H{z\Xr). 

Suppose that prior to the start of the process, the system has come to equilibrium with a thermal reservoir at 
temperature T and the reservoir has been removed. Therefore at t = 0 the microstate of the system will be treated 
as a random sample from a canonical distribution corresponding to the Hamiltonian PIa- From t = 0 to t = r, the 
system is described by a trajectory Zt evolving under Hamilton’s equations, and the work performed on the system is 
given by the difference between its initial and final energies, as per Eq. 2: 

W = Er-Eo = HB{Zr)-HAizo) . (4) 

For an ensemble of realizations of this process, the distribution of values of work performed on the system is 

pC^W) = JdErJ dEo P^(E,|Eo) P^iEo) <5(IV - P, + Eq) , (5) 

where Pa{Eo) is the probability distribution of initial energies, sampled from equilibrium, and P'^{Er\Eo) is the 
conditional probability distribution to end with a final energy Er, given an initial energy Eq. Let us consider these 
factors separately. 

P^ is simply the classical equilibrium energy distribution; 

PaIEo) = ^^e-P^°gA{EQ) , ( 6 ) 

where 

Zf = i y dz exp [-/3i7(z; A)] and gx{E) = ^ j dz5[E - H{z\X)] (7) 

are the classical partition function and density of states, /3 = {kBT)~^, and h is Planck’s constant. We have included 
the factors h~^ (which cancel in Eq. 6) for later convenience. Let us also define 

r2(£'. A) = J dz 9 [E — E[ (z; X)] = pdq , (8) 

which is the phase space volume enclosed by the energy shell E of the Hamiltonian i7, equivalently the integral of 
pdq around this energy shell. We then have 

1 dn 
hWE 


9x{E) 


(9) 
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(a) i = 0 



FIG. 1. An energy shell Ao of the initial Hamiltonian evolves under H{t) to the curve A = At- The curve B is an energy shell 
of the final Hamiltonian. 


The conditional probability distribution P'^{Et\Eq) is given by the expression 


P^{Et\Eo) = 


J dzp 6 [Eq - HAjzo)] S [Et - HBjzrjzo))] 
f dzo <5 [Eq - HaIzo)] 


( 10 ) 


where Zt{zo) denotes the final conditions of a trajectory that evolves from initial conditions zq. It is useful to imagine 
an ensemble of initial conditions sampled microcanonically from the energy shell Eq of the Hamiltonian EIa (see 
Fig. 1(a)). The swarm of trajectories that evolves from these initial conditions dehnes a time-dependent, closed curve 
At in phase space. At t = t, the points of intersection between At and the energy shell Et of the Hamiltonian PIb 
(black dots in Fig. 1(b)) represent the final conditions of those trajectories that end with energy Et, having begun 
with energy Eq. If we consider two nearby energy shells Et and Et -h dE, then P^{Et\Eq) dET is the fraction of 
trajectories whose final energies fall within this energy interval. 

Let us dehne the m ’th energy interval of the Hamiltonian P[a to be the range of energy values E satisfying 


mh < fl{E, A) < (m -|- \)h (11) 

where m = 0,1, 2 • • •. Furthermore, let us use a mid-point rule to assign a particular energy value E^ to each interval: 

= j> ^pdq= (^+^h . ( 12 ) 


Analogous definitions apply to the n’th energy interval of the Hamiltonian PIb, and the corresponding energy E^. 
From Eqs. 9 and 11 the width of the m’th energy interval of PI a is 


SEi 


For a smooth function of energy f{E), Eq. 13 leads to 


1 


9a{E^) 


dEofiEo) 


f{Et) 

9a{E^) 


(13) 


(14) 


where the integral is taken over the m’th interval. 

Readers will recognize Eq. 12 as the semiclassical quantization condition j> pdq = [m + (l/2)]/i, which provides an 
excellent approximation for the m’th eigenvalue of the quantum Hamiltonian Ha (Eq. 18 below) when m ^ 1. For 
convenience, we will later use the notation E^ to denote this eigenvalue, without making a distinction between the 
exact eigenvalue and its semiclassical approximation. 

The probability to obtain initial conditions zq within the m’th interval of Ha, when sampling from equilibrium, is 

Pi^(m)= / dF;oPf(Eo)«^e-^^- , (15) 

Jm 
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using Eq. 14. Similarly, the conditional probability to end in the n’th energy interval, given a representative initial 
energy in the m’th interval, is 





Ei) 


P^{E^\Et) 

9b{E^) 


(16) 


We will refer to this as the classical transition probability. The interpretation of P^(n|m) is straightforward. Imagine 
a swarm of trajectories evolving from initial conditions sampled from a microcanonical ensemble with energy 
(see Fig. 1(a)). At the end of the process, t = t, the fraction of these trajectories that fall into the energy window 
[E^,E^^,] is equal to P^{n\m) [66]. The classical work distribution can now be rewritten as a sum over initial and 
final energy intervals: 


P^iW)=Y^ / dE, dEoP^{Er\Eo)PS{Eo)S{W-Er+Eo) 

/ dA. / dEoP^iEr\Ei)P^iEo)S{W-E^+E^) 
= Y,P^in\m)P^{m)S{W-E^+E^) . 


(17) 


Finally, we point out that the approximations appearing in Eqs. 13 - 17 assume that the width of the m’th energy 
interval, 6Em, is very small on a classical scale. Formally, this assumption can be justified by taking the semiclassical 
limit, S —0, with all classical quantities held fixed. In this limit the density of states increases, as does the quantum 
number at a given energy. In practice the semiclassical limit is often simply identified with the condition m ^ 1. 


B. Quantum setup 


In the quantum version of this problem, F7(z; A) is replaced by the Hermitian operator 

A wavefunction tp{q,t) = evolves under the Schrodinger equation as H{Xt) is varied from H{Xo) = Ha to 

H{Xr) = Hb- The eigenstates and eigenvalues of Ha are specified as follows: 

= , HAlcj^i) = E^\cl>i) (19) 

with similar notation for Hb- Evolution in time is represented by the unitary operator Ut satisfying H{Xt)Ut = 
ihdUt/dt and Uq = I, where I is the identity operator. 

Using Eq. 1 and assuming thermal equilibration at t = 0, the work distribution is [15] 

pQ{W) = Y,P'^in\m)P^{m)6{W-E^+E^) . (20) 

n,m 


P^{m) is the probability of obtaining the m’th eigenstate of Ha when making the initial energy measurement: 

n 

Pa(.^) = , (21) 

with = ^m ^^P(~ldEm)- The quantum transition probability P^{n\m) is the conditional probability to obtain the 
n’th eigenstate of Hb upon making the final measurement, given the m’th eigenstate of Ha at the initial measurement: 


pQ{n\m)= U^\Ur\(l3t) = |(<(>f|V’r) 


( 22 ) 


The classical and quantum work distributions given by Eqs. 17 and 20 can now be compared directly. We note that 

Z^= f dEe-^^gA{E) « ^ 

m 


( 23 ) 
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hence P^{m) « P^{m) (see Eqs. 15 and 21). Thus what remains is to clarify the relationship between the classical 
and quantum transition probabilities, P^{n\m) and P^{n\m). In the following section, we compare these transition 
probabilities in a model system for which both the classical and quantum dynamics are simulated numerically. We will 
find that P^ and P^ are manifestly different (Fig. 2), but these differences mostly vanish after appropriate smoothing 
(Fig. 3). In Sections IV and V we develop a semiclassical theory to explain these features. 


III. NUMERICAL CASE STUDY: THE FORCED QUARTIC OSCILLATOR 

We begin with a model system, the forced quartic oscillator: 

Hiz;X) = ^+Xq^ . (24) 

We set M = 1/2 and h = h/2Tr = 1, and we vary the work parameter at a constant rate. A* = Aq +ut, from Aq = 1 = A 
to Xr = 5 = B, taking v = 50 and therefore r = 0.08. 

A number of groups have previously compared quantum and classical work distributions for a driven harmonic 
oscillator, for which analytical solutions are available. Specifically, Deffner et al [67, 68] have studied an oscillator 
with a time-varying stiffness; Talkner et al [69] have obtained the work distribution for an oscillator driven by time- 
dependent perturbations that are linear in q and p; Campisi [70] has studied changes in the Boltzmann entropy for 
forced quantum and classical oscillators; Ford et al [71] have investigated a harmonic oscillator with a cyclically driven 
stiffness; and Talkner et aZhave considered the sudden quench of a two-dimensional oscillator [26]. While the harmonic 
oscillator has the advantage of being analytically tractable, it is somewhat special in that its quantum dynamics can 
be reduced to its classical dynamics [72]. Here we have chosen the more “generic” quartic oscillator, which requires 
numerical simulations. 

To evaluate P^{n\m) we simulated 10^ Hamiltonian trajectories evolving under H{z;Xt). Initial conditions were 
sampled microcanonically, and final conditions were binned according to the energy intervals of Hb{z), as defined in 
Sec. HA. For our initial microcanonical ensemble, we chose m = 150, corresponding to Eq = = 1749.23. The 

results are shown by the dashed line in Fig. 2. The transition probability is a smooth function of n from nmin = 105 
to n„iax = 215, and zero outside this range. 

These sharp cutoffs reflect the fact that our initial conditions were sampled from a microcanonical distribution: 
nmin and rimax correspond to the minimal and maximal final energies that can be reached by trajectories launched 
with initial energies Eq = E^, as illustrated in Fig. 7 of Sec. V. The characteristic U-shape of the classical work 
distribution within the allowed region is related to the fact that the curves A and B shown in Fig. 1(b) become 
tangent to one another at A = Flmin and E = Ej^six (again, see Fig. 7). For an exactly solvable model in which the 
same U-shape appears, see Fig. 1 of Ref. [70]. 

pe(n|150),PTn|150) 



FIG. 2. Quantum (22) and classical (16) transition probabilities for the forced quartic oscillator. The solid black curve shows 
P‘^(n|m), while the dashed red curve shows P‘^(n|m), for m = 150. 
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For the quantum system, we expand a wavefunction evolving under H[Xt) 


/2M + \t(f as follows: 


= ^Cn(t) e 

n 



(25) 


Here H\4>n) = En\4>n) and the c„’s are expansion coefficients. The time-dependent Schrodinger equation ihip = H'4’ 
then produces the set of coupled ordinary differential equations, 


Cn 



dX 


= -A^ 

k^n 


{(j)n \q'^\ 4>k) 

Ek - En 


giC7Ti-7fc) 




(26) 


We integrated these equations numerically, from initial conditions c„(0) = 6mm to obtain 


P^{n\m) = |c„(t)P. 


(27) 


In Fig. 2 we plot the quantum transition probability as a function of the final quantum number n (solid line). 
Although a correspondence between P'^{n\m) and P^{n\m) is visually evident, the quantum and classical cases differ 
in two distinct ways: (i) the quantum probability oscillates rapidly with n, and (ii) P'^{n\m) shows tails that “tunnel” 
into the classically forbidden regions n < Umin and n > Umax- Both features trace their origin to the wave nature of 
the quantum system. We also see in Fig. 2 that P^{n\m) = 0 when n — m is an odd. This result follows from the 
fact that |g^| (f>k) = 0 when the parities of n and k differ, hence only transitions between states of the same parity 
occur under Eq. 26. 


^pe(k|150),^fTk|150) 

4=0 4=0 



FIG. 3. Accumulated transition probabilities for the forced quartic oscillator. The jagged black curve represents the quantum 
case, X)fe=o P'^{k\m), while the smooth red curve represents the classical case, P'^{k\m). 

In Fig. 3 we plot the accumulated transition probabilities X)fc=othereby smoothing 
out the rapid oscillations in the quantum transition probability. The close agreement observed in Fig. 3 suggests that 
Eq. 1 is indeed the appropriate quantum counterpart of Eq. 2, though distinctly non-classical features are visible in 
the interference and tunneling effects in Eig. 2. 

In Sections IV and V we investigate these issues analytically. We will find that both the agreement seen in Eig. 3 
and the quantum effects evident in Fig. 2 can be understood quantitatively, within a semiclassical interpretation in 
which quantum dynamics are approximated by classical trajectories bearing time-dependent phases. 


IV. SEMICLASSICAL THEORY 

In Sec. IV A we rewrite the classical transition probability P^{n\m) (Eq. 16) as a sum over trajectories that begin 
and end with energies and E ^, respectively (Eq. 39) . In Sec. IV B we use time-dependent WKB theory to derive 
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a semiclassical approximation {n\m) for the quantum transition probability P^{n\m). Our main result, Eq. 49, 
is expressed as a sum over the same trajectories that contribute to P‘^{n\m), only now each of these trajectories 
carries a phase, resulting in interference effects between different trajectories. When these interferences are ignored, 
we recover the classical transition probability (Eq. 50). When they are included, the interferences give rise to the 
oscillations observed in the quantum transition probability. In particular, interferences between symmetry-related 
trajectories account for the fastest oscillations observed in Fig. 2, as we discuss in Sec. IV C. 


A. Classical transition probabilities 


We begin with the conditional probability distribution P^{E\Eo) (Eq. 10), where we have dropped the subscript 
T from the final energy, for convenience. As discussed in Sec. II, we imagine a swarm of trajectories evolving in time 
under E[{z;Xt), with initial conditions sampled from the microcanonical phase space distribution 


. (.■F\= ^{Eo-HAiz)] ^ 6 [Eo-Ha{z)] 
Jdz6[Eo-'HAiz)] hgAiEo) 


(28) 


where Eq is a parameter of the distribution. The quantity P^{E\Eq) dE is the fraction of these trajectories that end 
with a final energy in the inhntesimal range {E, E + dE). 

As in Fig. 1, let At denote the time-dependent curve that evolves from the initial energy shell Ha = Eq, and let B 
denote an energy shell of the final Hamiltonian, Hb = E. The evolving curve At can be described by a multivalued, 
time-dependent momentum field where the index b labels the branches of At at the coordinate value q. E.g. 

in Fig. 1(b) the momentum field for Ar has two branches at q = qi. Similarly, the multivalued momentum held p^{q) 
describes the hxed energy shell Hb = E, whose two branches are 


p^iq) = ±^2M[E-V{q-B)] {b = ±). 


(29) 


For convenience, we will henceforth use the notation A (without a subscript) to indicate the surface At- The points of 
intersection between A and B, highlighted by dots in Fig. 1(b), represent the trajectories that contribute to P'"{E\Eq). 

Our swarm of trajectories at time t is described not only by the evolving surface At, but also by a probability 
density on that surface. If we project this density onto the g-axis, the projected density (g) is a sum over the 
branches of the surface At- 


/‘(g) = ^p^‘(g). 

b 


(30) 


The hxed microcanonical ensemble corresponding to Hb = E can similarly be projected onto a density 




1 

hgB{E) 


dp6[E - HB{q,p)] 


1 

hgaiE) 


E 


OHb ■ 
9P t 




(31) 


with OEb/Op = p/M evaluated at p = pf{q)- 

Let 1 = 1, 2, • • • AT be an index labeling the intersection points of A and B (e.g. AT = 4 in Fig. 1), and let bi denote 
the branch corresponding to the I'th intersection point: 

{qupi) = {quPbMi)) = {quPbMi)) ■ ( 32 ) 

With some abuse of notation, we have written bi rather than bf' and bf, though in general the Lth intersection point 
may occur at differently indexed branches of the two surfaces. 

To evaluate P^{E\Eq), we begin with the contribution Pf {E\Eq) from a single intersection point, 1. Fig. 4 depicts 
this intersection, together with an intersection point involving an inhnitesimally displaced energy shell Hb = E + SE. 
The highlighted line segment connecting these two points represents a set of trajectories with final energies between 
E and E + 6E, hence 

p^{qi)-\6q\=Pf{E\EQ)-\SE\. (33) 


Since we focus here on the contribution from a single intersection between A and B, we suppress the subscripts bi 
indicating the branches of these surfaces. By construction, the ratio 6E/Sq is the rate of change oi Hb with respect 
to g along the curve A: 






( 34 ) 
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FIG. 4. The filled circle indicates an intersection between the closed curves A and B, see Fig. 1. The open circle is an 
intersection point between A and a nearby energy shell Hb = E -\- SE. 


Since p®(g) specifies a surface of constant Hb, we have 

dp^ dp 
dq dq 


dHe/dq 
Hb dHs/dp 


hence 


d ^ . _4, dHe , dHBdp-^ dHB (dp^ dp’^ 


Combining Eqs. (31), (33), (34), and (36) we get 


Pf{E\E^)=p^{qi) 
Summing over all intersection points gives us 


6E 


Sq 


-1 


= hgB{E) p-^{qi)p^{qi) 


dp-^ dp^ 
dq dq 


K 


P^{E\E^) = hgB{E) J2p^M)pl{qi) ^ ^ 


-1 


If we now set ^ ? then Eq. 16 gives us 

K 


P^{n\m) = hJ2 P'^,iqi)PbMi) 


1^1 


dpf^ dpf^ 

dq dq 


-1 


(35) 


(36) 


(37) 


(38) 


(39) 


This is the main result of this subsection. 

We will use the term classically allowed region to denote the range ifmin "E E < ifmax within which the function 
P'^{E\E[)) does not vanish, representing all final energies that can be attained by trajectories with initial energy 
Eq. In the coarse-grained function P^(n\m), the classically allowed region is bracketed by the values n^ain and nmaxj 
denoting the energy intervals of Hb containing Emin and Emax- As illustrated in Fig. 7 below, the values Emin and 
Emax correspond, respectively, to the largest energy shell of Eg that fits entirely within the closed curve A, and the 
smallest energy shell oi Hb that surrounds the entire curve A. At these two energies, the curves A and B become 
tangent to one another, dpr^Jdq = dp’^Jdq, leading to divergences in Eq. (38), which are reflected in the two sharp 
peaks in P^{n\m) in Fig. 2. Similar divergences appear in a more familiar context, namely the density p^{q) that 
describes the projection of a fixed microcanonical ensemble Hb = E onto the coordinate axis (see Eq. (31)). In this 
case the divergences occur at classical turning points at which the velocity dHB/dp vanishes. 

In the example introduced in Sec. Ill, the reflection symmetry of the quartic potential implies that trajectories come 
in symmetry-related pairs: if {qt,Pt) is a solution of Hamilton’s equations under the time-dependent Hamiltonian 
H{q,p; At), then so is {—qt, —Pt)- This means that for every intersection point {qi,pi) between A and B, there will be 
another intersection point at {—qi,—pi), and the two will contribute equally to P'^{n\m). 
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B. Semiclassical transition probabilities 


Let us now evaluate P^{n\m) in the semiclassical limit. Appendix A provides a brief introduction to time-dependent 
WKB theory, and for further details we refer to the reviews by Delos [73] and Littlejohn [74]. Using Eq. A5, the 
wavefunction V’(<zU) = {Q\Ut\4'm) ^he n’th eigenstate of Hb can be written as 








^ oBf \ ■ B~ 


(40a) 

(40b) 


where the actions Sf* (q) and (q) generate the momentum fields 


.Aw X dSC' 


B( \ 


(41) 


and the integers and /r® are Maslov indices (see Appendix A). The densities and momentum fields in Eqs. 40 and 
41 are the same as those appearing in Sec. IV A. 

From Eq. (40) we get 


(^flV'T)=E / <^q\/Pb Pb' exp I - St,) - i {fJ-^ - Hb') 


b,b' 


(42) 


which appears as Eq. 37b in Ref. [66]. Using the stationary phase approximation to evaluate the integral, we obtain 
a sum of contributions from points qi satisfying the condition 




dq 


dq 


which is equivalent to 


P^iqi) =Pb'iqi) =Pi 


(43) 


(44) 


Hence {’Pnl’^’r) reduces to a sum of contributions arising from the K intersection points of the surfaces A and B, 
representing trajectories with initial and final energies E:^ and , just as in the classical case (Eq. 32). 

To determine the contribution from the Z’th intersection point, we perform a quadratic expansion around q = qi'. 


SbM) - St, (q) « AS'z -f ^Ki{q - qif 


with AS"/ = S^^{qi) - S^^iqi) and 


.. 5^4 4 d^S^ dp^ dpi 

- -^(qi) - -^(qi) - -^(qi) - ^im)- 


(45) 


(46) 


Performing the stationary phase integral then gives us: 

{(t>n\'^r) ~ E \lKPbi exp ij^^Si - iApi"^ 


1^1 

K 


dq exp 


IK-l , x2 


= H V Ph “P 


1^1 


127rh 


(47) 


,ZO'i 7 I -/4 _ 


= J2aiE 


where are evaluated at q = qi; Api = /i^ — /x®; and ai = sign(/Ci) = ±1. In Ref. [66], where similar calculations 
are performed, a quantity equivalent to our 7r(T//4 is identified as a Maslov-like phase (Ref. [66], p. 1037). 

Eq. (47) expresses (0 ^]'(/'t) as a sum over trajectories, each contributing an amplitude and a phase: 


ai = \ 2TTh 


Pbi Pbi 

l«^/l 


01 = - ^A/X; -p 


(48) 
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The semiclassical transition probability (n\m) « \ {4>n\ V’t)|^ is then given by 


P^(^{n\m) = 


K 

E 

1=1 


js. 


aie 


(49) 


The general form of Eq. 49 is not surprising. The path integral formulation of quantum mechanics tells us that solutions 
of the time-dependent Schrodinger equation can be represented in terms of sums over classical paths decorated by 
phases exp{iS/h), and when > 0 the dominant contributions to the sums come from those paths that are solutions 
of the classical equations of motion [75]. 

We expect the semiclassical transition probability function {n\m) to provide a good approximation to the 
quantum transition probability P'^{n\m), for large quantum numbers. For the forced quartic oscillator, Fig. 5 compares 
these functions, taking m = 150. P^ was computed as described in Sec. Ill, and P^^ was evaluated directly from 
classical simulations of trajectories evolving from an initial microcanonical ensemble. (The evolution of the multivalued 
action was obtained by integrating pdq — H dt along each trajectory, see e.g. Sec. 3.5 of Ref. [74].) We observe 
that the agreement is excellent, except in the vicinity of the boundaries between the classically allowed and forbidden 
regions. We will examine these boundaries in more detail in Sec. V. 

Now let us simplify Eq. 49 by simply ignoring the cross terms in the double sum. This is the so-called diagonal 
approximation [76-78], and it leads to 


{n\m) 


= 2TTnY,PhPh 


Z=1 


Z=1 


9Ph dpi 

dq dq 




(50) 


which is identical to the expression for the classical transition probability, ^*^( 711771 ) (Eq. 39). 


P‘2(n|150),p'^(n|150) 



FIG. 5. A comparison of quantum and semiclassical transition probabilities for our model system. The black solid curve shows 
P'^(n|m), while the red dashed curve shows P®‘^(n|m), as given by Eq. 49. All the parameters are the same as those in Fig. 2. 

The agreement between P^ and P^^ = \ ^ a/e*®' [^ seen in Fig. 5 suggests that the rapid oscillations of the quantum 
probability distribution can be understood in terms of phase interference between different trajectories. When these 
interferences are ignored, we recover the classical probability distribution, as we see in Eq. 50. Hence P^^ acts as a 
bridge to between the quantum and classical transition probabilities, 

P^{n\m) ~ P^^{n\m) P*" (n|m) , (51) 

and thus between the quantum and classical work distributions (Eqs. 20, 17). The first approximation in Eq. 51 
involves time-dependent WKB theory together with a stationary phase evaluation of integrals, and the second is the 
diagonal approximation in which interferences are ignored. 
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C. Interferences and symmetries 

Let us now examine the interferences in more detail. Recall from Sec. IV A that the trajectories contributing to 
Eq. 49 come in symmetry-related pairs. This allows us to write: 

K j 

(</>f IV'r) = XI = X + a;+je*®'+-^) (52) 

where J = i^/2 is an integer, and the Vth and {I + J)’th trajectories are related by symmetry: 

{qi+j,vi+j) = {-qu-vi) ■ (53) 

Here we have indexed the points {qi,pi) from I = 1 to K, according to the order in which they appear as we proceed 
clockwise around the energy shell B, first rightward along the upper branch and then leftward along the lower branch. 
Eq. 53 implies 


ai+j = ai , ai+j = ai . (54) 

Using the fact that the points I and I + J are located directly opposite one another on the closed curves A and B 
(Eq. 53), we obtain 

- Sf = ]^jj)Aq= [m+^TTh , Sf+j - Sf = j^pAq = {n ttK (55) 

and 

Ml+j — fJ-f = ^^f+J ~ mF = 1 I (56) 

hence 

ASi+j = ASi + {m- n)Trh , Api+j = Apn . (57) 

(In Eq. 55 we have invoked the semiclassical quantization condition, Eq. 12, together with Liouville’s theorem, which 
guarantees that p dq remains constant with time. In Eq. 56 we have used the fact that the Maslov index is 
incremented by +2 as one proceeds around the closed curve A or B, see e.g. Fig. 1 of Ref. [66].) Eqs. 48, 52, 54 and 
57 give us 




1^1 


1 + e' 


i{m—n)'7T 


(58) 


When m — n is even, the symmetry-related trajectories interfere constructively, but when m and n have opposite 
parities they interfere destructively and we get P^^{n\m) = 0. This provides a semiclassical explanation for the 
observation, noted in Sec. Ill, that the quantum transition probability vanishes when to — n is odd: the most rapid 
oscillations in P‘^(n|TO) arise from interference between symmetry-related trajectories. The same interpretation was 
given by Miller to explain angular momentum selection rules in the scattering matrix describing an atom impinging 
on a homonuclear diatomic molecule (see Sec. Ill C of Ref. [64]). 

The interference effects in , due to cross terms in Eq. 49, include contributions not only from pairs of trajectories 
that are related by symmetry (Eq. 53), but also from those that are not related by symmetry: 


P'^'^(n|TO) = X®?+ X + X , (59) 

1=1 l^k l^k 

where denotes a sum of symmetry-related pairs of trajectories, {qk,Pk) = —(qhPi), and is a sum over non- 
symmetry-related pairs, {qk,Pk) —{qhPi)- To illustrate these contributions separately, we constructed the quantity 
by omitting the sum on the right side of Eq. 59, and we similarly constructed P^^y^ by omitting The 
results are shown in Fig. 6. Pfy'^ displays a regular oscillation, as the symmetry-related trajectories are either exactly 
in phase or exactly out of phase, according to the parity of to — n. The oscillations in are less regular, as the 

phases 9i — 9k are not necessarily integer multiples of it. When all cross terms are included, both sets of oscillations 
combine to give the pattern in Fig. 5. For instance, the gap observed near n = 200 in Fig. 5 reflects the corresponding 
minimum in seen in Fig. 6(b). 
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Pf,‘i(n| 150),/■'■(nl 150) pSO^j„| 150).P^(n| 150) 




FIG. 6. Contributions to the interference pattern in {n\m). Figure (a) includes only interferences between symmetry- 
related trajectories, while (b) includes only interferences between non-symmetry-related trajectories. 


V. AIRY TAILS 

In Fig. 5 we see that that our semiclassical approximation fails around the boundaries of the classical allowed 
region at and Umax- Although the semiclassical transition probability derived in Sec. IV vanishes identically 
outside the classically allowed region (since only classical trajectories contribute to the quantum transition 

probability tunnels into the classically forbidden region, as seen in Fig. 5. In this section we construct a semiclassical 
approximation that describes this behavior. Our central results are given by Eqs. 64 and 69, and illustrated in Fig. 9. 

Similar tunneling behavior is observed in the wavefunctions of energy eigenstates, which exhibit tails that reach 
into classically forbidden regions [79], as well as in semiclassical treatments of the S matrix in molecular collisions 
[63]; in both cases the tails are approximated by Airy functions. The expressions that we obtain in this section will 
also be expressed in terms of Airy functions. 

The failure of P^'^{n\m) in the boundary regions originates in the fact that the curves A and B become tangent 
to one another at E = Ej^in and E = Ej^ax- In the example we have studied numerically the curves are tangent 
simultaneously at two points, due to the symmetry of the quartic potential, as shown in Fig. 7(b). However, it is 
useful first to discuss the generic case, where there exists only a single point of tangency, depicted in Fig. 7(a). For 
specificity we will discuss the upper boundary of the classically allowed region at Amax, but similar comments apply 
to the lower boundary at Emin- 

Let us consider the function 


ASiq,E)=S-^iq)-S^iq,E) 


(60) 


that appears inside the exponent in Eq. 42, and its derivative 


Ap{q,E)=p-^iq)-p'^{q,E) = ^ASiq,E) . (61) 

We have suppressed the subscripts b and b', and we have explicitly indicated that S'® and p® depend on the final energy 
shell E. In Sec. IV we evaluated the integral / dq exp(iAS/Ii) by performing a quadratic expansion of AS around 
a single intersection point qi between the curves A and B, where Ap vanishes (Eig. 4). As the energy E approaches 
Emax from below, two points of intersection between A and B coalesce at a point {qc,Pc), indicated in Fig. 7(a). Fig. 8 
illustrates the behavior of Ap(q, E) and AS{q, E) in the vicinity of qc and E-max- In particular, AS has two stationary 
points when E < Emax and none when E > Emax- The stationary phase integration performed in Sec. IV (Eqs. 43 
- 47) implicitly assumed well-isolated stationary points, but this assumption breaks down near E = Emax- In this 
situation we evaluate the integral / dg exp(iAS/Ii) by expanding AS to cubic rather than quadratic order. Leaving 
the details of the calculation to Appendix B, here we present the result: 


f 

i 

/2h'' 

1/3 

i 



/ dg exp 

-AS{q,E) 


1 exp 

J^^S{qn,E) 

Ai 

Emax) 


(62) 
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(a) (b) 

FIG. 7. Solid blue lines depict the surfaces A = At (as in Fig. 1). Dashed red lines are energy shells B of Hb- Diamonds and 
dots indicate the points of tangency between A and B at i7min and i?max, respectively, (a) In the absence of symmetries, we 
generically expect a single tangent point at Emin and i?max. (b) For the Hamiltonian that we studied numerically, Eq. 24, the 
symmetry of the quartic potential implies that both A and the energy shells B are symmetric under {q,p) —^ (—q, —p), thus 
points of tangency come in symmetry-related pairs. 


where 


k = 


Qq2 -^raayi) ; ^ \Qc^ -^raayi) 


(63) 


and Ai(-) denotes the Airy function. Combining this result with Eq. 42, we obtain the following expression, which is 
valid when « Emax'- 


P. 


SC,tail 


{n\m) = |(<^f|V'r)|^ =47rV^Pc 


2h\ 

W\) 


2/3 


Ai 


'k^/3hV3 




(64) 


where and the subscript 1 emphasizes that we have assumed only a single point of tangency between 

A and B. At the lower boundary of the classically allowed region, we obtain the same result but with Amax replaced 
by Amin in Eqs. 63 and 64. 

The Airy function decays rapidly for positive values of its argument, and oscillates for negative values: 




(65) 


for real C ^ 1- One can establish by inspection that v > Q > k when E > Amax and iz, fc > 0 when E < Amin, 
therefore the argument of the Airy function in js positive outside the classically allowed region. As a result 

the transition probability decays monotonically, penetrating the classically forbidden region over a characteristic skin 
depth \k^/3h^/3v-3\. 

Let us now consider the situation corresponding to the example of Sec. Ill, in which there are simultaneously 
two points of tangency between the curves A and A at A = Amax, and these two points are related by symmetry 
(Fig. 7(b)). Using the subscripts c and d to distinguish these two points, we have: 


(icPc) = (-%, -Pd) , Vc = -Vd = V 
Summing the contributions from these two points gives us 


kr = -kd = k 


A,B ^ A,B 
rc rd 


2h\ 


1/3 


2l/3l 


gi(ASc/fi-A/Jc7r/2) _|_ gi(ASd/?i-A/j,j7r/2) 


('/'n Ii/'t) - 27ryp^pf j Al 1 ^ 1 / 3 ^ 213 ^^ Amax) 

where ASc,d are evaluated at A = Amax- Following arguments similar to those leading to Eq. 58, we obtain 
gi(ASc//t—A/Xcir/Z) _j_ gi(ASd/fi—A/id'7r/2) _ ^i{ASc/h—/2) ^ _|_ ^i(m—n)'K 


( 66 ) 


(67) 


( 68 ) 
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Ap(q,E) AS(q,E)-AS(q„E) 




FIG. 8 . Behavior of Ap and AS for E ~ Smax- The open circles correspond to intersection points between A and B when 
E < Ema.x- These coalesce into a single point of tangency at E = iSmax (solid circle), and for E > i?niax, A and B do not 
intersect. 


Eqs. 67 and 68 lead to a result identical to Eq. 64, apart from a factor that captures the interference between the two 
parity-related trajectories: 


T-}SC,tciil / I \ T-)SC.tail / i \ /i 2 

F 2 {n\m) = ’ (n|m)-4cos 


(m — n)TT 


(69) 


with the subscript indicating two points of tangency between A and B. 

Fig. 9 compares the quantum transition probability with our semiclassical approximations given by (Eq. 49) 
and (Eq_ 59 )^ jn tbe central and boundary regions, respectively. The agreement is excellent, illustrating that 

the entire quantum work distribution can be understood in terms of contributions from individual classical trajectories, 
and the interferences between them. 


VI. DISCUSSION AND CONCLUSION 

The correspondence principle broadly states that classical mechanics is recovered from quantum mechanics in 
the limit of large quantum numbers. In textbook discussions this principle is often illustrated in terms of static 
properties, typically by comparing the coordinate space probability distribution of a harmonic oscillator eigenstate 
with the classical, microcanonical distribution at the same energy [79, 80]. Here, by contrast, we have considered 
the correspondence principle in a dynamic setting, comparing classical and quantum transition probabilities under a 
time-dependent Hamiltonian. In this setting, time-dependent WKB theory provides a set of tools for using classical 
trajectories to construct approximate solutions of quantum dynamics. 

We have applied these tools to study the relationship between the definitions of quantum and classical work given 
by Eqs. 1 and 2. Focusing on the transition probabilities P^{n\m) and P^{n\m) for systems with one degree of 
freedom, we have obtained three main conclusions, and have illustrated them with simulations of a driven quartic 
oscillator: 

• In the classically allowed region, the quantum transition probability can be approximated in terms of interfering 
classical trajectories (Eq. 49, Fig. 5). 

• When the interferences between these trajectories are ignored, the classical transition probability is obtained 
(Eq. 50, Fig. 3). 

• The tunneling of the quantum transition probability into the classically forbidden region is accurately described 
by a decaying Airy function, whose arguments are expressed in terms of classical quantities (Eqs. 64, 69, Fig. 9). 

The second conclusion, and Fig. 3 in particular, clarifies the sense in which classical mechanics is “recovered” in 
the limit of large quantum numbers (in this context), whereas the first and third conclusions describe inherently 
quantal effects - interference, tunneling - that can nevertheless be understood and approximated in terms of classical 
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p8(nll50),P^'^(n|150) 



FIG. 9. Comparison between the quantum transition probabilities evaluated by numerical integration of the Schrodinger 
equation (black, solid line), and the semiclassical approximations given by Eq. 49 in the central region 110 < n < 200 (red, 
dashed), and by Eq. 69 in the boundary regions 90 < n < 110 and 200 < n < 230 (green, dashed). 


trajectories. In view of similar analyses in atomic and molecular contexts (see e.g. Refs. [62-66]) these conclusions 
are not surprising. However, given recent interest in quantum nonequilibrium work relations, our results provide a 
timely investigation of the correspondence principle as it applies to the work performed on a system driven away 
from equilibrium. As we have shown, semiclassical mechanics provides a bridge between classical and quantal work 
distributions, and our conclusions provide some justification for the two-point measurement definition of quantum 
work (Eq. 1). 

It will be interesting to study whether our results generalize to systems with more than one degree of freedom. 
For Af-particle systems with ^ 1, we generically expect that the classical transition probability {n\m) will be 
bell-shaped rather than U-shaped, with Uniin and Umax found deep in the tails of P^{n\m). In this situation the 
quantum tunneling into the forbidden regions will be negligible, since both P^'{n\m) and P^{n\m) will be very small 
at Uniin and riniax- An obstacle to studying Af-particle systems analytically is the lack of semiclassical representations 
of energy eigenstates, analogous to Eq. 40b, for generic systems with multiple degrees of freedom. Progress might be 
made by studying two limiting cases: integrable and fully chaotic systems [81]. For systems whose classical dynamics 
are integrable, semiclassical expressions for energy eigenstates can be constructed using classical action-angle variables 
(see e.g. Eq. 11 of Ref. [82] or Eq. 3.5 of Ref. [83] or Eq. 37 of Ref. [84]). At the other extreme, for fully chaotic 
systems. Berry’s conjecture [85] suggests that energy eigenstates can be treated as Gaussian random functions of the 
configuration, q. In either case, the expression for the energy eigenstate might be combined with time-dependent 
WKB theory as a first step toward generalizing Eq. 42 and subsequent results to multi-dimensional systems. 
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Appendix A: Brief review of time-dependent WKB theory 


In time-dependent WKB theory, approximate solutions of the time-dependent Schrodinger equation are expressed 
in terms of classical constructions in phase space. In the brief summary that follows, we restrict ourselves to a single 
degree of freedom, and a Hamiltonian of the form H{q,p,t) = /2M -|- V{q,t). 

A wavefunction 'ip[q,t) is said to be in WKB form when written as the product of a slowly varying amplitude and 
a rapidly oscillating phase. 


ip{q,t) = A(g,t)exp 




(Al) 


or else as a sum of such terms, as in Eq. (A5) below. Locally, Eq. (Al) describes a wave train ijj oc exp(i/cg), with 
wave number k = h~^dS/dq. This implies a local momentum 

f)Q 

p{q,t) = hk = —{q,t) , (A2) 

and the Born rule implies a probability density 

p{q,t) = \A{q,t)\^ . (A3) 


These functions inherit their time-dependence from iplq^t)- A standard calculation, in which Eq. (Al) is substituted 
into the Schrodinger equation and h is treated as a small parameter [74], produces the equations of motion 





= 0 


dp 

dt 


dq 



= 0 


(A4a) 

(A4b) 


where terms of order h? have been neglected. Eq. A4a is the Hamilton-Jacobi equation. Eq. A4b is the continuity 
equation for evolution under a velocity held p{q,t)/M. 

To interpret these equations, let the function p{q, 0) describe a curve Ao in phase space, at an initial time t = 0. 
If the wavefunction ?/; is represented as a sum of terms of the form given by Eq. Al, then p(q,0) is multivalued, and 
the curve Ao has multiple branches, but for the moment we focus on a single branch. Ao is a Lagrangian manifold, 
determined by its generating function, S{q,0), via Eq. A2. Now imagine that this curve carries a probability density 
whose projection onto the coordinate axis is p{q,0). It is convenient to picture a swarm of particles, sprinkled along 
Aq so that the fraction of particles between q and q + dq is given by p{q, 0) dg. From these initial conditions, each 
particle in the swarm evolves under Hamilton’s equations, giving rise to a time-dependent curve At and density p{q, t). 
The Hamilton-Jacobi equation governs the evolution of At, through its generating function S{q, t), and the continuity 
equation governs the evolution of the projected density p{q,t). 

With these considerations in mind, approximate solutions of the time-dependent Schrodinger equation are written 
in the form 


■0(9,0 = X! \/^b(9A) exp 

b 


i 

h 


Sb{q,t) - ifib 


TT 

2 


(AS) 


The sum is taken over the branches of a Lagrangian manifold At whose generating function S{q,t) satisfies the 
Hamilton-Jacobi equation (Eq. A4a), and the projected probability density for each branch satishes the continuity 
equation (Eq. A4b). 

The quantities pb are Maslov indices. These determine the relative quantum phases of the various branches of At, 
and they differ from one another by integer values. (By convention, the overall phase of is adjusted so that the pbS 
themselves are integers.) See Ref. [74] for a more detailed discussion of Maslov indices, as well as a careful treatment 
of caustic points, where two branches of At meet and the manifold becomes “vertical” in the (g,p)-plane. 

At the semiclassical level of approximation, Eq. (A5) connects the evolution of a quantum wavefunction to that of 
a swarm of classical particles “surfing” on a Lagrangian manifold. This perspective provides a natural starting point 
for a comparison between quantum and classical work distributions. 
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As an important example of a WKB wavefunction, consider a manifold that is a single energy shell E oi a fixed 
Hamiltonian H{q,p), and consider a microcanonical distribution of initial conditions on this manifold. Under time 
evolution, neither the manifold nor the probability distribution changes, as each trajectory merely goes round the 
energy shell. The solution of Eq. A4a in this case is 

Sbiq,t) = Sbiq,0) - Et , (A6) 

hence by Eq. A5 the wavefunction simply acquires a phase Thus in the semiclassical limit, energy eigenstates 

correspond to microcanonical ensembles. Indeed, Eq. A5 in this case leads to the WKB approximation for energy 
eigenstates that is familiar from undergraduate textbooks on quantum mechanics, where it is typically derived in a 
different manner. Moreover, a proper treatment of the Maslov indices leads to the quantization condition j> p{q) dq = 
[m + {l/2)]h. 


Appendix B: Semiclassical transition probabilities in the boundary area 


Expanding the function AS{q,E) defined by Eq. 60 to cubic order, we get 

ASiq, E) « ASfiE) + AS'fiE) {q - qfi) + \aS'^{E) {q - q,f + ^ {q - q,f , (Bl) 

where ASfiE) = AS{qc, E) and primes indicate derivatives with respect to q. Note that here we expand AS around 
the coalescence point q = qc, where the surface A becomes tangent to the surface B (see Fig. 8). We now wish to 
evaluate the integral / dq exp{iAS/h), using Eq. Bl. 

Let us first consider the integral 


/ + 00 

dq exp 

-OO 




, f{q) = ao +aiq+^a 2 q'^ E^asq^ 


Rewriting f{q) as 


where 


c = —- 


Q!2 

Ctz 


f{q) = fio +filiq-c) +-fisiq-cf , 


_ aia2 _ a\ o _ ^3 

— ao -h —^ , Pi — ai — -— , P3 — ^ j 

03 6 ai 203 2 


we obtain 


r+cc 

I = / dy exp 


^ ( /3i2/+ 


1/3 

= f ^ ] Ai 


/3i 




using the integral representation of the Airy function. 


1 /■+°° 

Ai(C) = ^ y dt exp 




We will now use this result to evaluate 


l[E) = dq exp 


^AS{q,E) 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


for small H, applying the stationary phase approximation and treating E as a parameter of the integral. 
The cubic expansion for AS given by Eq. Bl gives us the coefficients 


ao{E) = AS{qc,E) 

(B8) 

ai{E) = AS'{qc,E) = Ap{qc,E) 

(B9) 

a2{E) = AS"{qc,E) = Ap'{q,,E) 

(BIO) 

a3iE) = AS"'iq,,E) = Ap"iq,,E) 

(Bll) 
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with Ap = p-^ — p'^ = AS' (Eq. 61). Note that 

ai{Ec) = a2{Ec) = 0 


(B12) 


since A and B are tangent at qc when E = Ec- The integral I{E) is now given by Eq. (B5), with /3o, /?i and 
obtained from Oq, ai, CX 2 and via Eq. (B4). This result simplifies if we consider how the /3’s behave in the vicinity 
oi E = Ec as h ^ 0. 

Let us define 


Pc=P'^{qc) =p^{qc,Ec) 

(B13) 


(B14) 

q2 A q2 B 

(B15) 

~ BE 

(B16) 


Now look at the argument of the Airy function in Eq. (B5), 


and consider a fixed range of C values, C- < 0 < C+i chosen so that the asymptotic approximations for the Airy 
function (Eq. 65) are accurate at both and Noting that C{Ec,fi) = 0 and expanding ({E,h) to first order in 
{E — Ec), we see that as /i —>■ 0, the fixed range [C-, C+] translates to a range of energies whose width scales like 
Thus we will be interested in energies that have the following scaling relation: 

e - ^2/3 , (B18) 


where e = E — Ec- To leading order in h, and expressing quantities in terms of e rather than E, we have 

C = ;^2/35i/3 > ’^'^ere a=^(e = 0) , b = Paie = 0) . (B19) 

Now recall that /?i = Oi — a^/^a^. In the vicinity of e = 0 (where ai = 02 = 0), we have Oi ~ e ^ and 
a^l2a^ ~ g 2 ^ hence the latter term can be ignored. We thus write 


therefore 


'^“1 /n\ 





Combining results and discarding terms that vanish as b, —>■ 0, Eq. (B5) finally gives us 




1/3 


exp 




Ai 




(B20) 


(B21) 


(B22) 


which is equivalent to Eq. 62. 
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